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Abstract 

Previous experimental work on a two-dimensional (2D) electron gas in a Si- 
on-sapphire device led to the conclusion that both conductivity and phonon 
drag thermopower S 9 are affected to the same relative extent by weak lo- 
calization. The present paper presents further experimental and theoretical 
results on these transport coefficients for two very low mobility 2D electron 
gases in 5— doped GaAs/Ga z Ali_ x As quantum wells. The experiments were 
carried out in the temperature range 3-7K where phonon drag dominates the 
thermopower and, contrary to the previous work, the changes observed in the 
thermopower due to weak localization were found to be an order of magni- 
tude less than those in the conductivity. A theoretical framework for phonon 
drag thermopower in 2D and 3D semiconductors is presented which accounts 
for this insensitivity of S 9 to weak localization. It also provides transparent 
physical explanations of many previous experimental and theoretical results. 

PACS: 73.20 Fz, 73.40 Kp 
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I. INTRODUCTION 



Weak localization (WL) refers to the quantum correction to the conductivity which arises 
because the phase of the electronic wavefunction is not randomized during elastic collisions. 
Through interference effects, this leads to a decrease in the longitudinal conductivity o xx 
which can be restored either by raising the temperature, thus increasing the frequency of 
phase-disrupting inelastic collisions, or by applying a low magnetic field which destroys the 
phase coherence of the electronic paths. WL has been thoroughly investigated^ over the last 
two decades but the effect on other transport coefficients has received much less attention. It 
was predicted! that the electrical conductivity a xx and the electronic thermal conductivity 
X xx would be related by the Wiedemann- Franz relation in the usual way for elastic electronic 
scattering, and the available experimental datJ on the longitudinal components in a quasi- 
2D system are in accord with this prediction. (We take the 2D electron gas (2DEG) to be 
in the xy plane and, unless otherwise stated, the magnetic field B along z). 

The interpretation of thermopower data is complicated by the fact that there are two 

«-» <-> 

contributions, diffusion S d and phonon drag S 9 . There is theoretical consensus that the lon- 
gitudinal diffusion thermopower S d x should be modified by WL,H the result being equiv- 
alent to the well-known Mott relation between S xx and a xx for elastic electronic scattering. 
However, in 2DEGs, which are the focus of the present paper, S 9 X is very large at low 
temperatures so that S xx is difficult to observe. This situation is quite different from that 
encountered in 3D systems. Introducing impurities into a 3D lattice to reduce the elastic 
mean free path l e of the electrons (so that WL is observable in a xx ) inevitably leads to a 
strong reduction of phonon drag because of the scattering of phonons by the same impuri- 
ties. In the 2D case, l e is reduced by placing impurities into or near the 2DEG, but under 
normal circumstances these impurities have no effect on the phonons in the substrate and 
so drag is not thereby reduced. 

As far as we aware, there is only one theoretical paper! concerned with WL and phonon 
drag. This dealt with S 9 X but gave a result only for the range of magnetic field B where the 
localization term in a xx is quadratic in B, which is appropriate only at very low B, and the 
situation at arbitrary fields remained unclear. 

Experimental data for a 2DEG in a Si-on-sapphire device were published by Syme et 
a/.i which indicated that S 9 X and affected equally by WL. The significance of this 

result is most conveniently discussed, as did the authors of that work, in terms of the 
phenomenological equation relating the electric current density J to the gradient of the 
electrochemical potential E' = E + V^i/e, where E is the electric field and fi the chemical 
potential, and the temperature gradient VT, i.e. 

J = aE'-eVT (1) 

+-* <— > 

where e is a thermoelectric tensor. The thermopower S is the experimental quantity which 
gives access to e and is defined by E = S VT (the chemical part V/z/e is not observed in 
experimental situations, e.g. see Gurevich et al. measured under the condition J = 0. 

The coefficient e has two additive contributions, diffusion e d and phonon drag e 9 , which are 

responsible for the two corresponding components S d and S 9 . 
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S§. x varies rapidly with temperature and it is difficult to observe deviations that might be 
due to localization. It is therefore necessary to measure S xx as a function of perpendicular 

magnetic field B. From Eq. (|l|) with VT y = we have S xx = p xx e xx + p yx e xy} where p = a . 
The very low mobility of the previous (and present) samples means that it is an excellent 
approximation to write S xx = e xx /a xx . If we denote the relatively small changes due to B 
as Aa xx , etc., then we have 

Agg, Ae xx Aa xx 

S xx € X x & xx 

and the data were consistent with the first term on the right hand side being zero. Indeed in 
a further publication^! it was assumed that this was precisely so and residual field dependent 
effects at the level of < 1% were analyzed in terms of WL corrections to S xx . 

For reasons which will be explained in detail in Section |TH this experimental result 
is unexpected. Basically we will show that e 9 xx is proportional to the electron-impurity 
scattering rate, l/r ei , whereas S xx is not. One can view WL to be a modification of r ei due 
to interference effects. Thus, we would have expected both e xx and to be affected by 
WL, but not S%. x . This is just the opposite of that found by Syme et alB 

In view of these considerations, we have reexamined the question both experimentally 
and theoretically. The present paper gives new experimental data on two samples of 2DEGs 
in 5-doped GaAs/Gai-^Al^As quantum wells. In contrast to the previous results, we find 
AS^ X ~ for B < 1 T, which implies Aa xx /a xx = Ae xx /e xx (to an accuracy of ~ 10%). 

For several years we have been using various physical pictures of phonon drag ther- 
mopower to help in understanding experimental data on 2D and 3D systems in a qualitative 
way, e.g. Refs. and [13|. In Section |T|, we outline a theory of phonon drag thermopower 
which puts these ideas on a quantitative basis. In particular, it provides arguments to sup- 
port the result that WL has no influence on phonon drag. However, we also show that the 
theory provides a useful framework to understand many previously published experimental 
and theoretical results in a transparent physical way. A brief, preliminary account of some 
of this work has been presented in conference form.0-3 



II. EXPERIMENTAL TECHNIQUES AND RESULTS 

Both our samples have been well characterized in previous work.@ The electron densities 
are n = 3.6 x 10 16 m -2 for sample 1 (SI) and n = 2.21 x 10 16 m -2 for sample 2 (S2). Using 
the 4K values of the conductivities, each has a mobility of p ~ 0.12 m 2 / Vs. Only a single 
subband is occupied in each sample. 

The thermoelectric properties of SI have been previously studied as a function of B over 
a wide temperature ranged At low temperatures, where S xx for this device is dominated by 
phonon drag, no variation of S xx with B up to about 1 T was seen to within the experimental 
accuracy of < 2% at low fields, even though changes ~ 5 — 6% of a occurred due to 
localization effects. 

In the present experiments a different cryostat was used with the sample mounting and 
associated instrumentation redesigned to increase the precision of the measurements by 
about an order of magnitude in the range 3-7K where phonon drag is dominant in S xx , and 
also to allow the direction of the magnetic field to be set either perpendicular or parallel to 



3 



the 2DEG. The low temperature limit was chosen to avoid any problems at low magnetic 
fields caused by the superconducting/normal transition of In which is used to provide thermal 
contact between the sample and cold sink. 

The experimental techniques are based on those published previously. @ All data were 
taken by dc methods. Very slow field sweeps were used for p xx , with a sample current for 
which self-heating was negligible. To measure the temperature and temperature gradient 
we used a matched pair of 22 kQ Philips surface mount devices which have very high sensi- 
tivity at these temperatures and no observable magnetoresistance (even up to 8T). During 
a measurement of S xx , B was held constant and the temperature gradient was repeatedly 
switched on and off by routing a fixed power either to the sample heater, which produced the 
temperature gradient, or to a matched heater mounted on the cold sink close to the sample. 
In this way the average power input to the sample stage remained constant and times to 
equilibrate were short. With temperature differences mostly in the range of ~100-300mK 
we had a resolution for VT X of < 0.2%. It is important to note that VT X was independently 
measured at each value of B. Thermoelectric voltages of 3.5-22 /xV were measured with 
a resolution of a few nV by an EM 11 nanovoltmeter. The voltages were reproducible to 
< 0.1% when VT X was switched on and off at fixed T. 

During the measurement of S xx the average sample T was found to increase slightly 
with B. The relative change was always largest at highest T; at 7.4K and 1.5T it reached 
a maximum of about 0.8% for SI, and 0.4% for S2. On lowering the field to 0.5T, the 
changes dropped to 0.2% for SI and 0.1% for S2. These changes are similar for both 
perpendicular and parallel B. Each data point has been corrected to bring them all to the 
same temperature. Interestingly, when this is done, the results are essentially the same as if 
we had used the zero-field value for VT X and the as-measured values for the thermoelectric 
voltages. This arises because the thermal conductivity of the substrate and S xx have similar 
T dependences. The estimated relative error of S xx as a function of B is typically ~ 0.3% 
and that of p xx is < 0.1%. There is a systematic overall uncertainty in S xx of about 10% due 
to possible errors in thermometer and voltage contact spacing, but this should not depend 
on T or B. 

The temperature dependence of S xx for the two samples is shown in Fig. 1 along with 
S xx which was measured for Sl@ but is estimated for S2. The increase in S xx for S2 
compared with SI is due to the lower carrier density and somewhat higher substrate thermal 
conductivity. The data for SI are in good agreement with those of previous measurements.^ 

The upper panels of Figs. 2 and 3 show data on the B dependence of AS XX /S XX = 
{S XX {B) - S xx {0))/S xx (0) and Aa xx /a 

XX \&xx 

(B) — cr xx (0)) / a xx (0) (taking a xx — l/p xx ) 
for the two samples. The plots are given in terms of In B because if WL is dominant we 
expect a xx oc In I? over a relatively wide field range. Although Aa xx /a xx is about 6% 
and 10% respectively for the two samples, the variation of AS XX /S XX is roughly an order of 
magnitude smaller over the same field range, though it is not zero within experimental error. 
These data are contrary to those previously publishedi where AS XX /S XX showed essentially 
the same behaviour as —Aa xx /a xx . 

Although there seems no doubt that WL is responsible for the behaviour of a xx , we 
briefly outline our analysis of the temperature and field dependence of these data to show 
that they do indeed behave as expected. We expect the Boltzmann value of a xx (which we 
write as na^k F l e where = e 2 /2ir 2 h, k F is the magnitude of the Fermi wave vector and 
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l e is the elastic mean free path) to be modified due to contributions from WL, Aa^. L , and 
electron-electron interaction, Aa xx , of the form0 

(Txx = n<r N k F l e + AaZ L + A(t^ = a N [7rk F l e -(1+g) Hh/Q] (3) 

where U = v F T{ is the inelastic mean free path given in terms of the phase coherence time 
Tj, and g is a measure of the strength of the electron-electron interaction. The temperature 
dependence arises from ij. If one assumes U oc 1/T P , then a xx = p(l + g)a^ lnT+ constant. 
Plots of a xx versus In T over the range 2-9 K for both samples gave good straight lines with 
slopes of 1.93 x KT 5 ^ 1 for SI and 1.85 x lO^lT 1 for S2. From the analysis of the B 
dependence given below we find p = 1.0; these results imply g = 0.57 for SI and 0.51 for S2, 
each ~ ±10%. 

Noting that Aa xx is 5-independent for low fields, the 5-dependence of o xx was fitted to 
the standard results of WL,0 i.e. 



&xx — &N 



nk F l e - gin + V> (~ + ^j) - * (\ + ^ 



(4) 



where ip is the digamma function. Using an expression for the conductivity in this form puts 
a strong constraint on the value for l e because it is required to fit simultaneously both the 
Boltzmann term and the correction terms. The usual practice has been to fit only the field 
dependent part in isolation. As shown in Figs. 2 and 3 this expression gives excellent fits to 
our experimental data in the range —0.05 T< B < 0.05 T, but at higher fields the fits become 
worse and are no longer acceptable for \B\ > 0.1 T. This agrees with expectations since the 
field dependence in Eq. (|j) is accurate only when 2eBl 2 e /Ti <C 1,0 which corresponds to 
\B\ <C 0.25 T in our case. Within experimental error = a/T and l e is precisely constant 
(using fits in the range — 0.05T< B < 0.05 T); a and l e are given in Table 1, along with 
kpl e . If electron-electron effects are ignored in the above (g = 0), the fitted parameters are 
changed only weakly. 

Our samples are in the dirty limit with h/r e ^> ksT, where ks is the Boltzmann constant 
and r e the impurity relaxation time. Under these conditions, Altshuler et a/Shave predicted 
that U should be given by 

where ep is the Fermi energy and r e is the electronic elastic relaxation time. This expression 
yields U = 51/T/zm for SI and 31/T/xm for S2, which are to be compared with the exper- 
imental values of 14/Tum and 21/T/im. The agreement is certainly not as good as that 
obtained by Choi et a/.H3 with GaAs heterojunctions, but the order of magnitude is correct, 
especially for S2. All these results confirm that our samples are dominated by WL, at least 
insofar as the B dependence of a xx is concerned. 

In the previous work! data were also taken with the magnetic field parallel to the 2DEG, 
where no change is expected nor observed for a xx , and no change was seen for S xx . We have 
also done this and typical data on S xx (with B || y) are shown in the lower panels of Figs. 2 
and 3 for both samples. Any changes are now at the level 0.2% which is within the range of 
the expected uncertainties. 
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III. THEORY 



In keeping with the experimental situation, we shall restrict our attention to the case of 
a 2DEG in which the electrons occupy a single subband and have an isotropic energy disper- 
sion. The generalization to multiple subbands is straightforward but will not be considered 
in this paper. Our treatment of the transport properties is based on the semiclassical Boltz- 
mann equation and is similar to the approaches used by Butcher's group@ and Lyo.iil It is 
useful to provide a brief but complete summary of the theory in order to clearly reveal the 
underlying approximations that are made. In addition, our intention is to bring out various 
features of the theory, not immediately evident in the previous theoretical work, which give 
useful physical insights into the general behaviour of the phonon drag component of the 
thermopower in degenerate semiconductors. A result of particular interest is the definition 
of an effective electric field which accounts for the effect of a nonequilibrium flux of phonons 
on the electrons. Once these semiclassical results have been established, we indicate the way 
in which the theory can be extended to include the effects of weak localization. 

In the presence of static electric and magnetic fields, the electron distribution function 
/k(r) satisfies the steady-state equation 



v k • V/ k - | (E + v k x B) • V k / k (>J k 



+ 0/k 



dt 

imp 



(6) 

ph 



where the terms on the right hand side account for impurity and phonon scattering. The 
impurity scattering term is given by 



Of 



= E IM 1 - A) - /k(i - /k<)] (7) 

imp k' 



where is an energy conserving impurity scattering rate. Similarly, the phonon scattering 
term is given 



dt 



= E LMi - /k)P k 'k - / k (i - /kOfldc] (8) 

ph k' 



with 



2vr 

p kk' = -E WMQ)? [N Q 5(ew - e k - hcu Q ) + (AL Q + l)S(e v -e k + hu Q )} . (9) 

n Q 

The variable Q represents both the three-dimensional phonon wavevector Q with compo- 
nents parallel (Qy) and perpendicular (Q±) to the 2DEG, and polarization index A; uiq is the 
phonon frequency and Nq is the nonequilibrium phonon occupation number. The term pro- 
portional to Nq corresponds to phonon absorption whereas the (ALq + 1) term corresponds 
to emission. The electron-phonon matrix element has the property M kk /(Q) = M£, k (—Q) 
and is proportional to <5 k ' 5 Q ||+k , due to momentum conservation in the plane of the 2DEG, 
and to the subband matrix element I(Q±) = J exp(iQ ^z)\4>(z)\ 2 dz. The detailed form of 
A^kk' (Q) depends on the nature of the interaction; the acoustic phonon deformation poten- 
tial and piezoelectric interactions are relevant to the GaAs/GaAlAs systems studied here. 
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Finally it should be noted that electronic screening of the electron-phonon interaction is in- 
cluded by dividing the bare electron-phonon matrix element by the dielectric function e(Qii) 
of the 2DEG. 

The solution to Eq. (§) is obtained by linearizing about a local equilibrium (le) Fermi 
distribution, i.e. 



/k(r) = /i e (T(r) )A i(r)) 



dey- 



$1 



(10) 



where T(r) = T + ST(r) is the local temperature and /x(r) = fi + 5/x(r) is the local chemical 
potential. The thermodynamic equilibrium distribution at temperature T and chemical 



potential ji will be denoted by /£ = (e /3 ( £k ~^ + l) -1 . 
is written in the form shown for later convenience, 
expanded as 



The deviation from local equilibrium 
Similarly, the phonon distribution is 



Nr 



I dN Q 



(11) 



where Nq = [eP^Q — l) -1 is the equilibrium Bose distribution at temperature T. We shall 
assume that the nonequilibrium phonon distribution Gq is established by a temperature 
gradient and is insensitive to the interactions of the phonons with the electronic system. 
However, phonon-drag thermopower arises directly from the interaction of the electrons 
with the nonequilibrium distribution of phonons. 

Substituting these forms of the electron and phonon distributions into Eq. (|6]), we obtain 
the linearized equation (details of the derivation of the phonon scattering terms can be found 
in the work of Cantrell and ButcheiS) 



<9/ k ° df k 



v k .VT^-ev k .E'^ + -v k xB.V k ^ k ^ 



+ 



imp 



dt 



(12) 



ph 



where 



dt 



imp 



k / C^k 



(13) 



and 



dt 



ph 



with 



(3 £ - $ k } [TMQ) + r k 'k(Q)] - P £ G (Q) FwiQ) 

k',Q k',Q 



r kk '(Q) = y |M kk ,(Q)| 2 / k °(l - / k ° K<K^ - e k - hue, 



r k 'k(Q)] (14) 



(15) 



The effective electric field E' in Eq. (12) is now the gradient of the electrochemical potential, 
E + V/i/e. The first sum on the right hand side of Eq. flI3|), which we call dh/dt\ p h t i, 
represents the scattering of the nonequilibrium distribution of electrons from the thermal 
distribution of phonons which has the effect of equilibrating the electrons to the lattice. The 
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second sum, dfk/dt\ p h,2, accounts for the effect of the nonequilibrium phonons in driving the 
electrons out of equilibrium. When considered in isolation, the first sum is responsible for 
the resistivity due to phonon scattering. However, for the samples of interest here, impurity 
scattering is the dominant scattering mechanism at low temperatures and this term can be 
neglected in comparison to Eq. ([13]). The equilibrium-phonon scattering term is nevertheless 
of interest as will be discussed later. 

We now concentrate on the nonequilibrium phonon scattering term, d f^/ dt\ p h,2i which 
is responsible for phonon-drag thermopower. It depends on the nonequilibrium phonon 
distribution which is obtained from the phonon kinetic equation 



v ph (Q) • VN Q 



dt 



coll 



N, 



(16) 



where the phonon relaxation time r v h is a parameter characterizing boundary and impurity 
scattering. In principle it is a function of Q but we will take it to be a simple Q-independent 
constant which can be obtained from the phonon thermal conductivity. With our earlier 
definition of Gq, Eq. ([II]) implies that 



G c 



T 



■TphYphiQ) ■ VT. 



(17) 



Within an isotropic Debye model, the phonon velocity v p h(Q) = s\Q, has distinct values for 
the longitudinal and transverse modes. 

The distribution function in Eq. (fT7| ) has the property G(Q) = —G(—Q), which is the 
asymmetry one would expect to see in the presence of a temperature gradient. Using this 



property, the nonequilibrium phonon scattering term in Eq. (|T4 ) can be written in the form 



dt 



= -/?£ G(Q) [r kk ,(Q) + r k , k (-Q)] 

ph,2 k',Q 



(18) 



and substituting Eq. fll7|) into Eq. fll8|) , we find 



dt 



ph,2 



k BT 2 £ Q 



E ■ Qh [r kk '(<3) + rv k (-Q)] 



(19) 



The appearance of the in-plane projection Qn in the sum assumes that the temperature 
gradient is in the plane of the 2DEG. 



To simplify Eq.(19) it is convenient to make use of the following identities 

ldf° 



f(e)[l-f(e + r^)]N°(u) 
/°(e-M[l-/WM 



(3 de 
ldf 

(3 de 



N°(lu) + f{e + huj) 
N°(uj) + 1 - /°(e - hu) 



(20) 



for the thermal factors appearing in the scattering rates r k k'(Q) defined in Eq. (|15|). With 
these results, together with the momentum conservation condition Qii = k' — k, Eq.(19) can 
be written in the form 



8 



dt 



«9/ k ° ™**aAa 



ph,2 



E 



v k 



VT 
~T~ 



(21) 



where = s\T p h is the phonon mean free path and t*(e) is an energy- dependent electron- 
phonon relaxation time defined as 



1 



x 



J2 ^i-^cob^ |M kk ,(g)| 2 

{ K + AO] ^ - ^ - + + 1 - /V)] ^' ~ e + • (22) 



2tt 



Here, # is the scattering angle between the initial wavevector k and the final wavevector k' 
and we have introduced the short-hand notation e = £ k and e' = e k '- Provided that Eq . (|i~7l) 
is an accurate representation of the nonequilibrum phonon distribution, Eq. fl2~T|) is a rigorous 
result for the phonon collision integral. 

If we now define the effective electric field 



E 



ph\ 



= E 



m*s\K\ VT 



T 



A eT ep( £ ) 



(23) 



Eq. ( pip can be written in the suggestive form 



0/k 



dt 



ev k • E pfc (e 



<9e k 



(24) 



In the context of the original Boltzmann equation in Eq. fll2|), the electron-phonon scattering 
term given by Eq. fl24|) (as stated earlier, we neglect the dfk/dt\ p h,i term) can be grouped 
together with the term in Eq. ( |T2"D arising from the actual electric field E'. In other words, 
the phonon-drag field, E p ^(e), proportional to the temperature gradient, has exactly the 
same effect on the nonequilibrium electron distribution as the actual electric field. Its 
dependence on energy through the electron-phonon relaxation time is only an apparent 
complication. Since the dominant impurity scattering mechanism is elastic, the impurity 
Boltzmann equation can be solved at a given energy e to define an energy-dependent impurity 
conductivity a(e; B). The contribution to the current density arising from the phonon-drag 
field can then be expressed in the form 



Iph 



de 



de 



er(e; B) E ph (e) . 



(25) 



It can be shown that this result for the current density is equivalent to that given by Zianni 
et a/.il and therefore leads to the same results for the thermoelectric tensor and phonon-drag 
thermopower. 

Our interest here is in the low temperature properties of a degenerate 2DEG. Since the 
impurity conductivity at low magnetic fields depends only weakly on energy, it can be taken 
out of the integral in Eq . (|25|) and we simply obtain the energy average of the effective field, 



>ph 



a{B) (E ph ) 



(26) 
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where cr(B) is the measurable conductivity appearing in Eq. (ffl). In view of Eqs. ( |2"3"D and 
(p6|), the thermoelectric tensor is 

?(B) = _ £ ^/j_W (27) 



eT \ t a 



where the average lifetime is given by 




Se ) r^e) ■ (28) 

The phonon drag thermopower tensor, defined as S 9 (B) = a (B) e 9 (B), is therefore given 
by S 9 (B) = S3 1 with 

We see that the phonon drag thermopower is a diagonal tensor independent of the magnetic 
field. This result is a direct consequence of the thermoelectric tensor being a scalar multiple 
of the conductivity. Equations (|23|) through fl29|) are the main results of the theory. 

An equation similar to Eq. ( p9f ) was first given by Herring!! for 3D non-degenerate semi- 
conductors, and is explained in terms of the following physical picture. The non-equilibrium 
phonons preferentially flow down the temperature gradient VT carrying a crystal momentum 
current oc AVT. Scattering of electrons by phonons leads to momentum being transferred to 
the electrons at a rate proportional to l/r ep , where r ep is the e-p momentum relaxation time. 
This acceleration of the electrons proceeds for a mean time determined by the e-i momentum 
relaxation time r ei , (we are assuming, as above, that r ep <C r ei ) at which point impurity scat- 
tering randomizes the momentum. The mean drift velocity of the electrons established in this 
way gives rise to the electric current contribution 3 p h- Thus J p h = -e 9 VT oc (r ei /r ep )AVT. 
Since the thermopower is measured with J = 0, a compensating drift current, J a = o~E, is 
established. With J p h + J a = and a oc r ei , the induced electric field is E oc AVT/r ep , and 
S 9 has the form given in Eq. (|29|) . 

Before proceeding, we would like to indicate the usefulness of the results we have obtained 
so far. In the Appendix we show that the electron-phonon relaxation time appearing in 
Eq. ( ppp is closely related to the electron-phonon transport lifetime that arises in the context 
of the phonon-limited mobility. The latter is defined by fi ep = e(l/r tr )^ 1 /m* where the 
transport lifetime is given in Eq. (A4) of the Appendix. At low temperatures, (1/t^) 



and (1/t^) are essentially the same, in which case the phonon-drag thermopower can be 
expressed as S 9 ~ —Y,x(sxA x /n^ p T) with = e(l/r^,) _1 /m*. This result has been applied 
to degenerate 2DEGs (e.g. see Ref. but it has been viewed as a semi-quantitative result, 
which indeed is the case for non-degenerate semiconductors. The fact that it is quantitatively 
accurate for the degenerate case (at least in the low temperature limit) was discovered 
empirically in recent experimental workEll concerned with the thermopower of 2DEGs and 
composite fermions, and was used there to evaluate \i ep for these systems. The importance 
of this connection is that the contribution to the resistivity from phonon scattering can 
be very small and difficult to measure at low temperatures, even in high mobility systems. 



10 



Through Eq. (|29|), the thermopower provides an alternate way to measure /i ep which can be 
much more accurate. 

Equation (p9|) shows that the longitudinal phonon drag thermopower S^ x should be 
independent of magnetic field, and that the phonon drag contribution to the Nernst- 
Ettingshausen coefficient should be zero. These results are contained in the theoretical 
work of Zianni et alB but were presented in a form which is not as physically transparent 
as that given here. The field-independence of S^. x has been demonstrated experimentally 



in a number of 2DEG systems, e.g. Refs. 16 and 25, though it is also evident in much of 



the earlier work. However in contrast to the theoretical prediction, it was found that S^ x is 
not zero for the same systems. Butcher and Tsaousidoucj have suggested anisotropy of the 
electrons and phonons as a possible reason for this discrepancy, but this explanation would 
also imply some field variation of Sjj. x . 

Although the above results were derived with 2DEGs in mind, it should be emphasized 
that they are equally valid in three dimensions when the same physical conditions prevail. 
In Eq. (^) for example, the summation over k' is simply interpreted as a three-dimensional 
sum, all electron-phonon matrix elements are evaluated using three-dimensional plane wave 
states and Qii is replaced by the total phonon wave vector Q. The subsequent analysis 
leading to the thermopower in Eq. ( |2T3| ) remains unchanged. It was recently demonstrated 
experimentally that S^ x = to high accuracy in a degenerate 3D semiconductor, even 
when S%, x is completely dominant in the longitudinal case,t3 and that S% x is independent of 
magnetic field. In this respect, the 3D situation conforms even more closely to the theoretical 
predictions. 

We now turn to the aspect of immediate concern in this paper, that of WL. The weak 
field dependence shown by our experimental results for the thermopower is consistent with 
the conclusions reached above on the basis of the semiclassical Boltzmann equation, even 
though the behaviour of the conductivity in our samples is manifestly nonclassical. We 
argue that these observations have a natural explanation in terms of a generalized Boltz- 
mann equation developed by Hershfield and Ambegaokar.0 These authors showed that the 
effects of coherent backscattering can be included with the addition of an extra term in the 
semiclassical Boltzmann equation. Although their final results are derived in the absence 
of a magnetic field, they claim that their approach can be extended to include a magnetic 
field and that all the standard WL results for the magnetoresistance can be reproduced. 
The significance of this is that WL effects are in principle accessible within an otherwise 
semiclassical approach. 

We accept the argument that a suitable modification of the impurity scattering term in 
Eq. (|5]) can be made which, in the absence of phonon scattering, leads to a conductivity with 
WL effects included, and now consider the additional effects of phonon scattering. This is 
one inelastic scattering mechanism contributing to the phase-coherence lifetime n appearing 
in the WL correction, though electron-electron scattering is usually the dominant mechanism 
at low temperatures. Of interest here are the additional effects of phonon collisions which 
impart momentum to the electrons and induce a nonequilibrium electron distribution. 

Syme et a/.0>§l argued that, since phonon scattering events are phase disrupting, and 
since phonon drag originates in these events, e 9 xx will not exhibit WL effects. We believe 
that this is incorrect. Within a semiclassical description of the electron dynamics, we have 
shown that the effect of phonon collisions is equivalent to an effective electric field in so far 
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as the subsequent current response is concerned. Whether the impulse is provided by an 
electric field or a phonon collision, the induced current is limited by impurity scattering and 
is proportional to the impurity conductivity in either case. The same conclusion should also 
apply in the WL regime when r ep 3> r ei . Once an impulse has been applied, an electron 
propagates in the presence of the impurities and the full effects of quantum interference are 
operative up to a time determined by t%. It is important to emphasize that the frequency of 
electron-phonon scattering events is unchanged when the phonon distribution is displaced 
from equilibrium by the temperature gradient. This means that the frequency of phase dis- 
rupting events is unchanged and the relevant conductivity determining the current response 
is the one including the WL effects. Thus with the use of the generalized Boltzmann equa- 
tion to incorporate WL effects, our final results for the phonon drag thermoelectric tensor 
given in Eq. fl2~T|) should still be valid with the simple replacement of the conductivity by 
the one including weak-localization corrections. 

As a final comment, we note that the above argument used for phonon drag can also 
be made in the case of the diffusion thermopower. The first term on the left hand side of 
Eq. ( |T2"D gives rise to the diffusion current 

J d = ~ e Jds (-^Pj (s - /i) a(e; B) ^ . (30) 

Using a conductivity with WL corrections in this expression leads to the accepted form of 
the diffusion thermopower which, in contrast to the phonon-drag thermopower, does exhibit 
WL corrections.@~@ 

To summarize, we have argued that the role of phonon drag in establishing an electric 
current is equivalent to an applied electric field in the impurity-dominated regime, with 
the consequence that the phonon drag thermoelectric tensor is proportional to the impurity 
conductivity. The WL corrections within the conductivity also appear in the thermoelectric 
tensor, resulting in a cancellation of WL effects in phonon drag thermopower. These results 
are essentially in accord with the present experimental data. 



IV. DISCUSSION 

The origin of the discrepancies between the present and previous work with regard to S!j. x 
in a perpendicular magnetic field is not clear. The fact that the previous data were obtained 
on a MOSFET and the present data are for GaAs/Gai-^Al^As quantum wells should be 
irrelevant if WL is the cause of the field dependence of S%. x . This in itself is significant 
because if the differences are real, it implies that something other than weak localization is 
involved in at least one of these experiments. One possibility might be that the previous 
samples had kpl e ~ 2 — 4, though WL theory is expected to be valid only for k F l e 3> 1. 
The present samples have kpl e > 10 (Table 1). Much less is known, both experimentally 
and theoretically, about the situation at small k F l e where the system is approaching strong 
localization. Ashe et al^ have probed a xx in this limit but were unable to find satisfactory 
agreement with WL, even with extensions to the theory. How phonon drag might behave in 
this limit is unknown. 

There are also some features of the experimental technique used in the previous work that 
might be cause for concern. Syme et a/.B0§l used an ac method for which the temperature 
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difference could be measured only when B = and it was assumed to remain unchanged 
at finite B. Also, although constantan potential leads were used to minimize heat losses 
when thermal conductivity was measured, these were replaced by Cu leads for S xx . The 
authors argued that this gives isothermal conditions in the direction transverse to the ap- 
plied temperature gradient. However, the use of Cu leads implies a significant heat leak 
to the surroundings and unknown thermal boundary conditions on the 2DEG; this might 
be particularly important with an AC measurement where the heat flow is continuously 
changing. 

Although the change we see in S^. x is much smaller than that in a xx , nevertheless it is 
outside experimental error, especially for S2, and we must ask how it arises. An obvious 
possibility is that it is due to diffusion. The effect of WL on S xx , say AS XX , is believed to 

bM 

AS d xx „ Hk/l e ) 



St. nk F l e 



(31) 



For our samples at 3K (7K), \n(k/l e ) /nk F l e = 0.08 (0.06) for SI, and 0.12 (0.01) for S2. 
From Fig. 1 we see that for SI (S2) the measurements show that S^. x /S xx varies between 
about 25-70 (40-65) for the temperature range 3-7K. This suggests that changes due to this 
cause would be ~ 0.3% at 3K and ~ 0.1% at 7K. These estimates agree with the observations 
for SI, where AS XX /S XX is very small, but seem rather too low to account for the variations 
for S2. We would expect the effect to decrease as temperature increases thereby decreasing 
the relative importance of S xx , but the experimental results are not unambiguous on this 
point. We also note that these effects might be much larger, and basically unknown, for the 
samples of Syme et a/.§0 with kpl e ~ 2 — 4. 

Another possibility is that we may be observing the effects of anisotropy in the electron 
and phonon system. As we noted in Section [ITI], Butcher and Tsaousidou^ have suggested 
this as the origin of the finite S^ x observed in 2DEGs. If this is correct, it will be accompanied 
by a field dependence of S%. x , and this should be essentially independent of temperature, not 
inconsistent with the observations. 



V. CONCLUSIONS 

Our measurements show that the longitudinal phonon drag thermopower S!j. x is essen- 
tially independent of magnetic field, which implies that it is independent of weak localization, 
a result contrary to earlier work. We have given a theory of phonon drag based on the semi- 
classical Boltzmann approach which shows that Sfj. x is independent of magnetic field. Using 
the results of Hershfield and Ambegaokar0, we have argued that this result remains valid in 
the presence of weak localization. Our theory also gives results which are useful in under- 
standing the origin and behaviour of phonon drag thermopower in 2D and 3D degenerate 
systems. 
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APPENDIX A: 



It is apparent that Eq. ([19]) bears a strong similarity to the first phonon scattering term in 
Eq. (|T4| ) which is dependent on the nonequilibrium electron distribution. This similarity can 
be exploited to provide a relation between the phonon drag driving term given by Eq. (|19"D 
and the electron-phonon momentum relaxation time which determines the phonon-limited 
mobility. This latter quantity is obtained by solving the following transport equation 



ev k ■ E 



<9e k 



-p Y, i $ k' - $ k} [TMQ) + rVk(-Q)] 

k',Q 



(Al) 



which describes the transport of electrons scattering exclusively from phonons, that is, in 
the absence of impurities. It should also be noted that we have not included a magnetic field 
or thermal gradient in this equation. This of course is not the physical situation of interest 
in Eq. (|12|), but rather represents an auxiliary problem which identifies a useful transport 
property that is relevant to the phonon-drag calculation. In addition, a discussion of Eq. (|A1|) 
will allow us to make contact with earlier work&H on the calculation of electron-phonon 
mobilities. 

The right hand side of Eq. ( |Al"| ) follows from Eq. ( |i"4"|) with the replacement Q — > —Q in 
the r k / k (Q) term. This puts it in the same form as the corresponding term in Eq. (|19|). To 
proceed with the solution of (|A1|) , we make use of the identities in Eq. (|20|) . The appearance 
of the Fermi function derivatives allows the transport equation to be simplified as 



ev k • E 



2n V- 



$k}|M kk ,(Q)| 2 { 



N° Q + /V)l 8{e? 



hue 



+ 



[N Q + l-f(e')]S(s'-s + hu Q )}. (A2) 



We now look for a solution to this equation which has the form $ k = — ev k • ET tr {s). We 
shall refer to r tr {e) as the electron-phonon transport lifetime to distinguish it from the 
electron-phonon relaxation time defined in Eq. (22). With this ansatz, Eq. (|A2|) becomes 



2 71 " r k' cos 6* 

-r2^{ T tr{e) ; T tr 



k',Q 



k 



(e')}\M kk ,(Q)\ 2 { 



N°n + f(e> 



5(e' 



Tiojf 



+ 



[N Q + l~f°(e')]6(e'-e + hu Q )}. (A3) 



This is an integral equation for T tr (e) and is equivalent to Eq. (22) in Ref. [ |32 
Eq. (|A]|) is obviously very similar to Eq. 



The form of 



q). It simplifies if we now make the quasielastic 
approximationH which assumes that the energies of interest are near the Fermi energy ef 
and Huq <C £p. We can then replace e' by e in the r-dependent factor in Eq. (|A3|), and we 
obtain the following explicit expression for the relaxation time: 



Ttr(e) 



- cos i 

k'Q 



|M kk ,(Q)p 



X 



{[N° Q + f(e>)]6(e' 



e - hun) + 



[iV°+l-/V)] 5(e'-e + hu Q )} 



(A4) 



It can be shown that this expression for T tr (e) is identical to Eq. (39) in Ref. [ |32| where the 
quasielastic approximation is also invoked. The summation over Q includes the polarization 
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index A so that l/r tr (e) = Y^x^-/ T ir( £ )- ^ * s d ear that if the quasielastic approximation 
is also made in Eq. @), l/r^ (e) and l/r^ T {e) would in fact be identical. We expect any 
quantitative differences to be minimal at low tempertures where the two relaxation times 
can be used interchangeably. 

Finally, we can make contact with the average electron-phonon scattering rate which was 
derived by Stormer et a/.EH on the basis of the work of PriceH. If Eq. ( |A4|) is averaged over 
energy with respect to the weight function (—df°/de), we obtain the result 

— ) = x Et 1 - cos^|M kk ,(Q)| 2 /3{/ (e')[l " / V + ^q)]N° q 

Ttr' n k / g 

+ /°(s')[l-/V-^Q)]W + l)}. (A5) 



In arriving at this expression, we have again used the identities in Eq. ( PUD as well as 
the identity f°(e - huj)[l - f°(e)\N (iu) = /°(e)[l - f°(e - huj)](N (uj) + 1). The result 
actually given by Stormer et a/Eil follows from Eq.( |A5[ ) by performing the Qy sum with 
the understanding that k = k', so that Q» = 2/c' sin(#/2). The phonon-limited mobility 
defined as \i ep = e(l / T tr )~ l / m* was used by Stormer et to interpret their observation 
of a transition into the Bloch-Griineisen regime of the mobility of a 2DEG. Eq. (|A5[) is 
applicable at arbitrary temperatures, but a simpler expression can be obtained in the limit 
of low temperatures. Since the derivative of the Fermi function is then sharply peaked at the 
chemical potential /x, the average in Eq. (|A5|) may be replaced by the value at the chemical 
potential if T tr (e) is a slowly varying function of energy. In this case we have (see Eq. (|S3j)) 

^ = f E(i - ^ K^pn^) {^ k+Q " - " - hujQ) + 6{£ ^ - " + • (A6) 

In obtaining this result, we have made use of the momentum conservation condition 
|^kk'(Q)| 2 = l^(Q)| 2 ^k',k+Qn to perform the sum over k'. The magnitude of k in Eq. (|A6|) 
is now restricted to the Fermi wavevector. A similar approximation in Eq. (|28|) would allow 
the average lifetime defined there to be replaced by Eq. 
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TABLES 



TABLE I. The various parameters for the samples. l e and U = a/T were obtained from the fits 
of the data to the field dependence of a xx . 





n(10 16 m~ 2 ) 


H (m 2 /Vs) 


l e (nm) 


a (/imK) 




SI 


3.6 


0.134 


35 


14 


20 


S2 


2.21 


0.144 


42 


21 


13 
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FIGURES 



FIG. 1. The open circles give the measured thermopower S xx of the two samples. The lines are 
the diffusion components S xx . In the case of SI, S!~ x = —0.41 /jW/K is a measured quantity taken 
from Ref. 13. For S2, S xx = — 0.67 //V/K is estimated from that of SI assuming S xx oc 1/n. 

FIG. 2. A selection of the experimental data on the field dependence of a xx (o) and S xx (•) at 
various temperatures for SI. The results are averages of ±B data, though any dependence on the 
direction of B is weak. The upper and lower panels give data with B perpendicular and parallel 
to the 2DEG respectively. For clarity the data have been offset vertically by multiples of 3%. The 
solid line through the data for a xx in the upper panel is a theoretical fit according to Eq. (||). 



FIG. 3. The same as Fig. 2 except these data are for S2. 
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